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[ 1 ] Severe smog episodes over China in January 2013 
received worldwide attention. This air pollution was dis- 
tinguished by heavy loadings of fine particulate matter 
and SO 2 . To characterize these episodes, we employed the 
Ozone Mapping and Profiler Suite, Nadir Mapper (OMPS 
NM), an ultraviolet (UV) spectrometer flying on the Suomi 
National Polar-orbiting Partnership (SNPP) spacecraft since 
October 2011. We developed an advanced algorithm to 
quantify SO 2 in the lower troposphere and achieved high- 
quality retrievals from OMPS NM, which are charac- 
terized by high precision, ~0.2 Dobson Units (DU; 1 
DU = 2.69xl0 16 molecules/cm 2 ) for instantaneous field 
of view SO 2 data and low biases (within ±0.2 DU). 
Here we report SO 2 retrievals and UV aerosol index 
data for these pollution events. The SO 2 columns and 
the areas covered by high pollutant concentrations are 
quantified; the results reveal for the first time the full 
extent (an area of ~10 6 km 2 containing up to 60 kt of 
SO 2 ) of these episodes. Citation: Yang, K., R. R. Dickerson, 
S. A. Carn, C. Ge, and J. Wang (2013), First observations of 
SO 2 from the satellite Suomi NPP OMPS: Widespread air pol- 
lution events over China, Geophys. Res. Lett., 40, 4951-4962, 
doi: 10. 1002/grl. 50952. 

1. Introduction 

[ 2 ] The recent winter air pollution episodes over China, 
as publicized by news reports of sustained dense smog that 
engulfed Beijing in January 2013, are among the worst in 
history, rivaling the deadly Great Smog of London in 1952 
[ Thorsheim , 2009] and the suffocating photochemical smog 
of Los Angeles [Haagen-Smit, 1952] in the 1940s to 1970s. 
The largest sources of air pollutants in China are coal com- 
bustion and motor vehicle exhausts [He et ah, 2002], which, 
respectively, were also the major culprits of the historical 
London and LA smogs. China’s fast growing economy, as 
evidenced by its rapid industrialization and urbanization, is 
predominantly powered by its abundant supplies of coal, 
which accounts for about 70% of China’s total energy con- 
sumption [ U. .S'. Energy Information Administration, 2012]. 
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In addition to greenhouse gas (carbon dioxide, CO 2 ) emis- 
sion, coal burning releases significant amounts of air pol- 
lutants, including sulfur dioxide (SO 2 ), nitrogen oxides 
(NOx = NO + NO 2 ), carbon monoxide (CO), mercury (Hg), 
and fine particles (i.e., smoke). The gaseous pollutants, SO 2 
and NOx, are usually oxidized in the air and combine with 
water vapor to form sulfuric and nitric acid aerosols, which 
may be further mixed with the emitted smoke, ammonia, as 
well as dust from the deserts of Mongolia and China; this 
type of mixture is usually the main component of winter 
smog over eastern China [Tao et al., 2012]. 

[ 3 ] Anthropogenic S0 2 is the key signature of air pollu- 
tion in China, and its column density serves as a proxy for 
the concentration of pollutants coemitted from fossil fuel 
combustion. This is especially true in the winter months, 
when relatively low-tech heating by burning coal contributes 
a large proportion of air pollutants. These pollutants are usu- 
ally released into the planetary boundary layer (PBL) but 
may be subsequently lofted into the free troposphere by the 
winter cold fronts from the Siberian high [Wang et al, 2004; 
Dickerson et al, 2007] and transported over long distances 
[Hsu et al., 2012]. A synoptic perspective on the severity 
and extent of air pollution can be obtained using spacebome 
ultraviolet (UV) instruments, capable of detecting its sig- 
nature component, SO 2 , in the lower atmosphere [Eisinger 
and Burrows, 1998; Krotkov et al, 2006; Carn et al, 2007; 
Lee et al, 2008; Nowlan et al, 2011] and providing daily 
global coverage to monitor its evolution and transport. In 
this paper we report the first S0 2 measurements using the 
new UV Ozone Mapping and Profiler Suite (OMPS) Nadir 
Mapper (NM) instrument, launched into orbit on 28 October 
2011 onboard the Suomi National Polar-orbiting Partnership 
(SNPP) spacecraft. We describe the first observations of SO 2 
from SNPP/OMPS, quantify the daily loadings, estimate the 
spatial extent, and characterize the evolution of the January 
2013 air pollution events over China. 

2. Retrievals From OMPS NM 

[ 4 ] Suomi NPP/OMPS NM is a nadir- vie wing hyper- 
spectral instrument that measures backscattered ultraviolet 
radiance spectra in the 300-380 nm wavelength range at a 
resolution of 1 nm (defined as the full width at half max- 
imum (FWHM) of the spectral response function) and a 
sampling rate of 2.4 pixels per FWHM. It is in a Sun- 
synchronous, low Earth orbit, with a local ascending (north- 
bound) equator crossing time at 1:30 P.M. Contiguous 
daily global coverage is delivered using a two-dimensional 
charge-coupled device (CCD) detector array that covers a 
2800 km cross-track swath (110° field of view, FOV) with 
a nadir pixel size of 50 x 50 km, attributed to the her- 
itage of the Total Ozone Mapping Spectrometer missions 
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Figure 1 . Maps of SO 2 vertical column (VCD) density and aerosol index (AI) from OMPS iterative spectral fitting (ISF) 
retrievals for 3 days in January 2013. (a-c) The SO 2 loading, maximum SO 2 VCD, and its corresponding location displayed 
above each SCL map. The loading is computed as the total SO 2 mass observed within the dashed square area on the map, 
with lower left comer at 101.9°E, 21.8°N and the upper right comer at 125. 0°E and 44.9°N. (g-i) The NCEP reanalysis 
850 hPa geopotential height (GPH in kilometers) at 06:00 UTC for these 3 days. Note that the full swaths of the OMPS 
observations are used in producing these daily maps, without filtering based on surface albedo, cloud fraction, or viewing 
angles. Only a few pixels with low retrieval quality are excluded from the maps. 


[Heath et al, 1975; Fleig et al, 1982]. The improved 
performance of OMPS NM includes high signal-to-noise 
ratio (e.g., ~1500 at 310 nm) radiance measurements using 
a single detector array to cover the entire measured spectral 
range; hence, it is uniquely suitable for SO 2 measurements. 
In contrast, the predecessors of OMPS NM, including the 
Global Ozone Monitoring Experiment (GOME) [Burrows 
et al., 1999] on the European Remote Sensing Satellite- 
2 (ERS-2), the Scanning Imaging Absorption Spectrometer 
for Atmospheric Chartography (SC1AMACHY) on Envi- 
ronmental Satellite [Bovensmann et al., 1999], the Ozone 
Monitoring Instmment (OMI) on Aura [Levelt et al, 2006], 
and GOME2 [Munro et al., 2006] on MetOp-A and B, 
measured earthshine radiance using three or more spectral 
channels to cover a broader wavelength range at higher 
sampling frequencies and narrower spectral resolutions. For 
SCIAMACHY, OMI, and GOME-2, there is a channel split 
near 310 nm. Consequently, the measured radiance spectra 
of these instruments in the channel-overlap region (extend- 
ing a few nanometers either side of 310 nm) are degraded 
and are in general not used in SO 2 retrievals, even though 
this region is highly sensitive to SO 2 in the PBL. We also 
stress that OMPS NM is currently the only operational UV 
satellite instmment providing contiguous daily global cover- 
age at low- to midlatitudes, due to data gaps in the OMI and 
GOME-2 measurements. 


[5] To take advantage of the hyperspectral UV measure- 
ments and reduce the impact of measurement noise, in order 
to extract more reliable geophysical information, we refine 
and apply the iterative spectral fitting (ISF) algorithm devel- 
oped by Yang et al. [2009, 2010] to OMPS NM radiance 
spectra between 308 and 345 nm. This spectral range is 
divided into two subwindows in the iterative fitting process, 
308-333 nm and 333-345 nm. The former is used to derive 
the amounts of O3 and SO 2 constrained by the respective 
climatological O3 [Bhartia and Wellemeyer, 2002] and SO 2 
[Lee et al, 2009] profiles, while the latter is used for deter- 
mination of spectral surface reflectivity or cloud fraction. 
Specifically, we use the Mixed Lambert-Equivalent Reflec- 
tivity (MLER) model [Ahmad et al., 2004] to account for 
the effects of partial cloud and/or aerosol detected in the 
OMPS instantaneous FOV (IFOV), and the MLER cloud 
pressure is simultaneously retrieved by adjusting it until 
the calculated Rotational Raman Scattering filling-in factor 
[Joiner et al., 1995] matches the corresponding measure- 
ment. This pressure is usually lower than that of the ground 
surface and represents the level of an opaque surface that 
reproduces the combined reflection of the ground surface 
and scattering particles (clouds and/or aerosols) in the atmo- 
sphere. The lower the MLER pressure, the higher the altitude 
of these scattering particles in the atmosphere. Using this 
simultaneously derived cloud pressure and cloud fraction 
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facilitates the proper treatment of the partial shielding of 
trace gasses located below the scattering particles, as well as 
the enhancement of trace gasses located at higher altitudes. 

[6] The precision of the PBL SCL data from the OMPS 
ISF retrieval is about 0.2 Dobson Units (DU; 1 DU = 
2.69 xlO 16 molecules/cm 2 ), estimated from the standard 
deviation of background SCL measurements. However, it is 
more difficult to quantify the ISF SO 2 retrieval accuracy, 
which is mainly determined by two factors: the vertical dis- 
tribution of SO 2 and the presence of cloud and/or aerosol. A 
mismatch between the climatological and the true SO 2 pro- 
files, as well as a difference in the modeled and the actual 
mean paths of photons backscattered to the satellite sensor, 
can lead to retrieval errors, especially in the presence of 
cloud and/or aerosols. The ISF retrieval using OMI data has 
been validated against in situ aircraft measurements in China 
[He et al., 2012]. It was found that ISF retrievals yield accu- 
rate quantification of SO 2 when the concentration is high and 
is more accurate than retrievals using the band residual dif- 
ference method [Krotkov et al., 2006]. Frequently, the mass 
loading of a geographic region is aggregated from many 
satellite IFOVs. In this case, because the errors associated 
from profile shape assumptions and cloud/aerosol effects 
may cancel out in the summation, the uncertainty in the total 
is usually much smaller than that of the individual measure- 
ments and is estimated to be less than 10% when the total 
SO 2 mass loading is greater than 10 kilotons (kt). 

[7] In addition to O3 and SO 2 columns, the ISF retrieval 
provides useful information on atmospheric conditions and 
the underlying surface at the time of observation. The UV 
aerosol index (Al), proportional to the slope of spectral 
reflectivity [Yang et al., 2009 and 2010] in an IFOV, is also 
obtained from the retrieved MLER parameters. A positive 
Al value indicates the presence of UV absorbing aerosols: 
dust, smoke, or both. Thus, UV satellite measurements can 
potentially provide coincident information on O3, SO 2 , NO 2 , 
and particulate matter, which represent four of the six criteria 
air pollutants identified by the US Environmental Protection 
Agency. 

3. Results: SO 2 and Smog Episodes 

[8] The wealth of information retrieved from OMPS NM 
measurements and its daily global coverage provide a unique 
synoptic perspective on episodic air pollution events. Here 
we demonstrate this capability of monitoring global atmo- 
spheric environments with the first quantification of air 
pollution over China in January 2013. 

[9] The air quality related geophysical parameters are 
displayed in Figure 1, which show maps of SO 2 verti- 
cal column density (VCD) in DU and Al over East Asia, 
retrieved from OMPS observations for 3 days in January 
2013. The SO 2 maps reveal the spatial distribution of this 
air pollutant, one of the main components of smog. Simi- 
larly, the Al maps (Figures Id— If) display the geographic 
distribution of UV-absorbing aerosols. These patterns sug- 
gest that both local aerosols (smoke and urban/industrial 
aerosols) and dust blown from the Taklimakan and Gobi 
deserts of western China contribute to the aerosol load in 
eastern China. The 3 days selected for Figure 1 represent 
dates with heavy (Figure la), light (Figure lb), and moderate 
(Figure lc) SO 2 loadings. For the day with a high load- 
ing, the map (Figure la) shows the widespread distribution 


of SO 2 , blanketing a large portion of eastern China; a high- 
pressure system (Figure lg) that suppressed atmospheric 
dispersion in the region likely contributed to heavy smog. 
In contrast, on the days with medium and low loadings, hot 
spots of air pollution remain, such as in the Sichuan Basin 
and Beijing area (surrounded by mountains to the west, 
north, and northeast). These hot spots occurred persistently 
during this month, indicating the strength of local sources 
and/or effects of local topography and meteorology. How- 
ever, under the synoptic conditions that are dominated by 
the midlatitude trough (Figures lh and li), large-scale high 
loadings of SO 2 are unlikely to occur. The maps also reveal 
the typical trajectories of air pollution outflow: the north- 
east flow path (Figure lb) via high midlatitude large eddy 
systems and the east flow path (Figure lc) via the prevalent 
subtropical-midlatitude westerly winds. 

[10] To chart the evolution of air pollution in January 
2013, we show in Figure 2 the daily S0 2 loadings 
(Figure 2a), the extent of heavily SO 2 polluted areas 
(Figure 2b, full bars), and smog-covered areas (Figure 2b, 
sub (red) bars). The smog-covered areas are subsets of 
heavily SO 2 (>0.5 DU) polluted regions, in which the UV 
Al is also high (>1). The time series in Figure 2 shows four 
distinctive, synoptic-scale air pollution episodes. The first 
one was the longest, lasting for 12 days from the begin- 
ning of January. During the first episode, high SO 2 loadings, 
often mixed with high levels of UV-absorbing aerosols, cov- 
ered a large portion of eastern China. Additional large-scale 
episodes occurred in 14 to 16 January, 21 to 24 January, and 
from 28 to 30 January, respectively. There was a more local- 
ized smog episode, affecting the Beijing and Tianjin regions 
and their vicinity, which occurred in 18 and 19 January. 

[ 11 ] Daily satellite observations provide snapshots of SO 2 
spatial distributions, not only influenced by the source of 
emissions but also strongly affected by the local meteoro- 
logical conditions and topography. The recurrent satellite 
overpasses allow time averaging of these instantaneous mea- 
surements, yielding a stable spatial distribution that relates 
more closely to the locations and fluxes of the emission 
sources [Fioletov et al, 2011]. Figure 3a shows the average 
S0 2 VCD map derived from OMPS NM observations with 
less than 50% cloud cover (i.e., MLER cloud fraction <0.5) 
in the IFOVs over East Asia for January 2013. The high- 
est average OMPS VCDs (>1 DU) are all located just south 
and southeast of Beijing, distributed within three provinces 
(Shandong, Hebei, and Henan), indicating that S0 2 emis- 
sions from these provinces are the major contributors to air 
pollution in January. These high concentrations of SO 2 were 
likely intensified by synoptic anomalies, including a weaker 
Siberian high and a localized high-pressure anomaly over 
northeast China (Figure 3c), which resulted in persistent 
higher relative humidity, higher temperature, and weaker 
winds that favor the formation of smog. Moderately elevated 
SO 2 concentrations are also evident in Liaoning, Jiangsu, 
Shanxi, and Sichuan/Chongqing provinces, signifying the 
locations of additional sources of air pollutants. Lastly, 
isolated and low-level enhancements of SO 2 above the back- 
ground are detected in Hubei province, Ningxia Region, and 
Urumqi in Xinjiang Region, representing weaker emissions 
sources that contributed to the widespread distribution of 
these air pollution episodes. 

[12] To evaluate the OMPS SO 2 retrievals, we show 
in Figure 3b the corresponding average OMI S0 2 map, 
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Figure 2. (a) Observed SO 2 masses (kilotons) and (b) SO 2 covered areas (million square kilometers) within the square- 
bounding box shown in Figure 1, for each day (except Sunday) in January 2013. Results on Sundays are not shown, because 
OMPS NM is operated in a high-spatial resolution mode, with a reduced spectral coverage. The daily mass is the sum of all 
SO 2 values in the bounding box for the day, but areas are calculated for high SO 2 (>0.5 DU) regions, i.e., the sum of the 
yellow and red zones in the corresponding map in Figure 2c. (c) The red zones are defined as smog-covered areas, in which 
both S0 2 and the UV aerosol index (AI) are high (i.e., S0 2 > 0.5 DU and AI >1), and the yellow zones are high S0 2 and 
low AI areas (SO 2 > 0.5 DU and AI < 1). 
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Figure 3. (a) The average S0 2 VCD map measured by 
SNPP/OMPS NM with MLER cloud fraction less than 
50% over East Asia for January 2013, (b) the correspond- 
ing S0 2 map composited with the lower tropospheric S0 2 
column (ColumnAmountS02_TRL) of the operational OMI 
S0 2 product, and (c) the 500 hPa geopotential height (GPH) 
anomaly and 700 hPa wind vector anomaly (with respect 
to the 1980-2010 means) for January 2013. Letters “L” 
and “H”, respectively, indicate the negative and positive 
geopotential anomalies. The highest OMPS S0 2 VCD is 
1.65 DU, located at (118.3°E, 36.8°N), and the OMI S0 2 
VCD is 0.8 DU at this location. Note that the OMPS map 
includes observations from its full swath, while the OMI 
map (Figure 3b) is produced with cross-track samples 1 to 
21 (out of 1 to 60) to avoid the effect of the row anomaly 
(see http://www.knmi.nl/omi/research/product/rowanomaly- 
background.php). 

derived from the lower tropospheric (TRL) S0 2 VCDs 
(ColumnAmountS02_TRL in the operational OMI prod- 
uct, OMS02). The OMI TRL S0 2 VCDs are derived using 
the linear fit (LF) algorithm [Yang et al., 2007], with a 
near uniform a priori S0 2 profile from the ground to about 
5.5 km above sea level. Similar spatial distributions of ele- 
vated S0 2 are exhibited in Figures 3a and 3b, with the 
OMI TRL S0 2 VCDs mostly lower than the correspond- 
ing OMPS ISF values. The lower TRL S0 2 columns are 
expected due to the generally higher center of mass of the 


a priori profiles relative to the climatological profiles [Lee 
et al., 2009] used in OMPS ISF retrievals. This comparison 
also shows that OMI detects more of the S0 2 transport away 
from the sources, as illustrated in Figure 3b with elevated 
S0 2 values distributed further east, away from the coast of 
China, primarily due to OMI’s smaller pixel sizes that are 
more capable of seeing through broken clouds. These con- 
sistent spatial distributions of S0 2 , obtained using different 
(ISF and LF) algorithms from independent satellite measure- 
ments, provide an indirect validation of the new OMPS NM 
retrievals. 

4. Summary 

[13] We have shown that OMPS NM provides high- 
quality PBL S0 2 measurements, characterized by low biases 
(within ±0.2 DU) in both along- and cross-track direc- 
tions of the orbital swaths and by high precision (~0.2 DU) 
for its IFOVs. The OMPS measurements are about 7 times 
more precise than the current operational OMI PBL S0 2 
retrievals, which have a ~1.5 DU noise level [Krotkov 
et al, 2008]. The improved S0 2 measurements are made 
possible by the unique instrument design that provides high- 
quality spectral coverage of the 300-380 nm range using a 
single CCD detector array and by the stable instrument per- 
fonnance. Both the design and performance facilitate the 
accurate characterization and effective stray-light correc- 
tions of its hyperspectral radiance measurements, which in 
turn allow the ISF algorithm to take advantage of the high 
S0 2 sensitivity region below 310 nm. The broader consistent 
high-quality spectral coverage, coupled with advances in 
retrieval algorithm and radiance correction techniques, now 
permit daily global observation of S0 2 , with an improved 
detection limit (precision of ~0.2 DU) for S0 2 in the PBL 
for its IFOVs. The OMPS NM measurements can detect 
and quantify much smaller S0 2 amounts (~0.01 DU) with 
enhanced spatial resolution when averaging the repeated 
samplings of the same region over time. These improved 
measurements allow more accurate quantification of atmo- 
spheric S0 2 from volcanic emissions and anthropogenic pol- 
lution to provide new information on the relative importance 
of these sources for climate studies. 

[14] Using the high-quality OMPS data, we were able 
to quantify synoptic-scale air pollution events over China 
in January 2013. The cumulative S0 2 mass measured by 
OMPS is 1065 kt (=1.065xl0 12 g = 1.065 Tg). This rep- 
resents the lower limit of total S0 2 loading within the 
bounding box (23.1° x 23.1° dashed box in Figure 1) for 
this month, as S0 2 obscured from the satellite FOV by thick 
clouds and emissions on Sundays (on which OMPS operates 
in a mode with reduced spectral coverage) are not measured. 
This observed S0 2 mass is about half of the total S0 2 emis- 
sion (~2 Tg) [Yang et al., 2010] from the 2008 Kasatochi 
eruption (among the largest in the past decade), highlight- 
ing that coal combustion in China continues releasing large 
amounts of S0 2 into the atmosphere, and contributes greatly 
to the country’s air pollution today. 

[15] The OMPS NM is designed primarily for measur- 
ing O3 in the upper atmosphere. As demonstrated here, 
this well-characterized hyperspectral instrument also pro- 
vides sensitive observations of trace gasses in the PBL. We 
are developing new capabilities to extract greater and more 
accurate information from these observations, so that the 
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quality of OMPS trace gas products is sufficient to match or 
surpass those of the Aura/OMI environmental data records 
(EDR). We are on track to achieve this goal, as the next 
OMPS NM instruments, with an enhanced spatial resolution 
(12 x 12 km at nadir), will fly on the NOAA and NASA Joint 
Polar Satellite System (JPSS) spacecraft: JPSS-1 (scheduled 
to launch in early 2017) and JPSS-2. 
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